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ABSTRACT 

We  present  the  first  complete  far-infrared  spectrum  (43-197  pm)  of  M82,  the  brightest  infrared  galaxy 
in  the  sky,  taken  with  the  Long  Wavelength  Spectrometer  of  the  Infrared  Space  Observatory  (ISO).  We 
detected  seven  fine  structure  emission  lines,  [O  i]  63  and  145  pm,  [O  m]  52  and  88  pm  ,  [N  n]  122  pm, 
[N  m]  57  pm,  and  [C  n]  158  pm,  and  fitted  their  ratios  to  a  combination  starburst  and  photo¬ 
dissociation  region  (PDR)  model.  The  best  fit  is  obtained  with  H  ii  regions  with  n  =  250  cm-3,  an  ion¬ 
ization  parameter  of  10-3,5,  and  PDRs  with  n  =  103,3  cm-3  and  a  far-ultraviolet  flux  of  G0  =  102,8.  We 
applied  both  continuous  and  instantaneous  starburst  models,  with  our  best  fit  being  a  3-5  Myr  old 
instantaneous  burst  model  with  a  100  Me  cutoff.  We  also  detected  the  ground-state  rotational  line  of 
OH  in  absorption  at  119.4  pm.  No  excited  level  OH  transitions  are  apparent,  indicating  that  the  OH  is 
almost  entirely  in  its  ground  state  with  a  column  density  ~4  x  1014  cm-2.  The  spectral  energy  distribu¬ 
tion  over  the  long-wavelength  spectrometer  wavelength  range  is  well  fitted  with  a  48  K  dust  temperature 
and  an  optical  depth,  tDus1  oc  2_1. 

Subject  headings :  dust,  extinction  —  galaxies :  individual  (M82)  —  galaxies :  ISM  — 
galaxies :  starburst  —  infrared :  galaxies 


1.  INTRODUCTION 

Infrared  luminous  galaxies  emit  an  infrared  luminosity 
comparable  to  or  greater  than  their  optical  luminosity.  The 
starburst  galaxy  is  an  infrared  luminous  galaxy  converting 
its  molecular  interstellar  medium  (ISM)  into  stars  at  a  rate 
that  can  not  be  sustained  for  a  Hubble  time.  The  radiation 
from  these  new  stars  is  reprocessed  into  infrared  radiation 
by  the  dust  in  their  parental  molecular  clouds.  Other  pos¬ 
sible  energy  sources  for  infrared  luminous  galaxies  are 
active  galactic  nuclei  or  the  energy  of  shocks  resulting  from 
galaxy  interactions.  Recent  studies  (e.g.,  Genzel  et  al.  1998) 
support  the  hypothesis  that  most  of  the  brightest  infrared 
luminous  galaxies  are  predominantly  powered  by  recently 
formed  massive  stars;  however,  the  simultaneous  presence 
of  active  galactic  nuclei  and  active  star  formation  in  some 
galaxies  shows  that  both  processes  can  occur  in  the  same 
phase  of  the  evolution  of  luminous  infrared  galaxies. 

As  a  result  of  its  proximity  (3.63  Mpc;  Freedman  et  al. 
1994)  and  moderate  infrared  luminosity,  M82  is  the  bright¬ 
est  galaxy  in  the  infrared.  Because  of  its  recent  episode  of 
star  formation  (e.g.,  Rieke  et  al.  1993),  it  is  considered  to  be 
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the  prototypical  starburst  galaxy.  Sometime  around  108 9  yr 
ago,  M82  experienced  a  close  encounter  with  M81  (Yun, 
Ho,  &  Lo  1993).  This  gravitational  interaction  likely  pro¬ 
duced  the  bar  seen  (Telesco  et  al.  1991),  which  provided  the 
large  streaming  motions  necessary  to  funnel  large  amounts 
of  molecular  gas  to  the  galaxy’s  center.  This  large  reservoir 
of  molecular  gas,  ~2x  108  M0  (Wild  et  al.  1992),  provides 
the  fuel  for  the  ongoing  starburst. 

Studies  of  the  M82  nucleus  in  the  mid-infrared  (Telesco  et 
al.  1991)  saw  strong  off-center  hot  spots.  Using  the  velocity 
profiles  of  the  M82  far-infrared  fine  structure  lines,  Lord  et 
al.  (1996)  modeled  two  hot  spots  as  well.  These  hot  spots 
may  be  star  formation  sites  resulting  from  cloud-cloud  colli¬ 
sions  that  occur  because  of  orbital  crowding  near  the  inner 
Lindblad  resonances  (Kenney  et  al.  1992).  Alternatively, 
Satyapal  et  al.  (1997)  modeled  the  nucleus  of  M82  as  a 
region  of  outward  star  propagation,  with  the  hot  spots 
being  the  most  recent  and  brightest  areas  of  star  formation. 

Relatively  insensitive  to  extinction,  far-infrared  spectros¬ 
copy  can  provide  a  unique  probe  of  infrared-bright,  dust- 
obscured  galaxies  like  M82.  Line  ratios  may  be  used  to 
constrain  physical  parameters  and  tend  to  be  less  vulner¬ 
able  to  calibration  uncertainties  as  well.  It  is  also  in  the 
far-infrared  that  some  of  the  most  important  cooling  lines  of 
the  ISM  of  galaxies,  the  [O  i]  63  pm  and  [C  n]  158  pm  lines, 
appear.  A  close  examination  of  M82  will  provide  a  template 
for  future  comparisons  to  possible  starburst  galaxies, 
including  those  at  high  redshift. 

Our  observations  and  the  analysis  used  on  the  M82  far- 
infrared  spectrum  are  described  in  §  2,  and  §  3  presents  our 
fine  structure  line  results.  We  discuss  line  ratio  model 
fitting,  OH  absorption,  and  fitting  of  the  entire  spectrum  in 
§  4;  §  5  presents  our  final  conclusions. 

2.  OBSERVATIONS  AND  ANALYSIS 

We  present  Infrared  Space  Observatory  (ISO)  long- 
wavelength  spectrometer  (LWS;  Kessler  et  al.  1996;  Clegg 
et  al.  1996)  grating  mode  (43-196.7  pm,  X/AX  ~  200)  obser¬ 
vations  of  M82.  These  observations  were  taken  as  part  of 
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the  Long  Wavelength  Spectrometer  Consortium’s  Extra- 
galactic  Guaranteed-Time  Program.  All  observations  were 
made  in  1996  May  on  two  separate  orbits.  The  central  posi¬ 
tion  was  R.A.  (2000)  =  9h55m5283  and  decl.  (2000)  = 
69°40'45'.'9.  Our  main  observations  consisted  of  20  full 
grating  scans,  with  a  total  integration  over  all  scans  of  8  s  at 
each  wavelength  position  spaced  at  one-quarter  of  the 
resolution  element.  The  spectral  resolution  was  0.29  /mi  for 
the  43-93  n m  range  and  0.6  /mi  for  the  80-196  /mi  range. 
Also,  in  the  final  averaged  spectrum,  the  spectral  resolution 
between  80  and  93  /im  lies  between  these  two  values, 
resulting  from  our  averaging  of  the  overlap  of  two  detectors 
that  are  detecting  different  orders  of  the  grating.  We  also 
took  a  spectrum  offset  by  Aa  =  +  0!8  and  A<5  =  +17'  from 
the  center  of  the  galaxy  in  order  to  measure  the  back¬ 
ground.  This  offset  observation  had  the  same  integration 
time  but  was  over  a  more  limited  wavelength  range.  All  data 
used  in  this  paper  were  processed  through  the  LWS  Pipe¬ 
line  Version  7. 

The  spectra  were  flux-calibrated  with  respect  to  Uranus 
(Swinyard  et  al.  1998).  The  individual  detector  scans  were 
calibrated  to  within  +4%-6%  of  each  other,  based  on 
overlapping  detectors,  with  the  notable  exception  of  detec¬ 
tor  LW4,  covering  wavelengths  150-170  /mi,  about  15% 
lower.  Several  reliable  detectors  were  chosen,  and  the  rest  of 
the  detectors  were  multiplicatively  shifted  so  that  overlap¬ 
ping  points  would  match.  The  LWS  beam  is  roughly  80" 
FWHM,  but  it  does  show  some  variation  with  detector, 
ranging  from  65"-85"  FWHM  (Swinyard  et  al.  1998). 

We  performed  the  postpipeline  analysis  with  the  ISO 
Spectral  Analysis  Package  (ISAP).  Even  after  pipeline  pro¬ 
cessing  the  spectrum  remains  contaminated  with  bad  data 
point  “  glitches  ”  from  cosmic  ray  hits,  which  we  cleaned  by 
plotting  several  scans  on  top  of  each  other  and  looking  for 
the  characteristic  fast  rising  and  slowly  falling  glitch  shape. 
All  bad  data  points  were  then  removed  as  well  as  any  suspi¬ 
cious  and  widely  deviating  points  following  the  glitch, 
where  the  detector  showed  memory  effects.  Roughly  20%  of 
the  data  were  discarded  as  unusable. 

Raw  data  also  contain  a  fringing  pattern  that  is  believed 
to  arise  from  interference  between  the  beam  coming  off  the 
field  mirror  and  another  beam  reflected  by  the  substrate 
holding  that  mirror.  Fringing  is  only  a  minor  effect  for  the 
shorter  wavelengths  ( <  70  /mi) :  the  fringing  is  weaker,  and 
the  spectrum  is  inherently  noisier.  At  wavelengths  longer 
than  100  /mi,  fringes  approach  ~  5%  of  the  continuum  and 
endanger  our  ability  to  measure  accurately  the  continuum 
level  as  well  as  to  identify  possible  weak  emission  and 
absorption  lines.  The  fringe  is  sinusoidal  in  wavenumber 


with  a  well-measured  period  of  3.54  cm-1.  We  used  the 
sine-wave-fitting  algorithm  for  extended  sources  in  ISAP  to 
divide  out  the  fringes.  The  defringing  has  little  effect  on  the 
majority  of  line  fluxes,  altering  the  weak  145  and  122  /mi 
lines  by  only  5%.  Only  the  OH  absorption  at  119  /mi 
changes  significantly,  becoming  33%  weaker.  Finally,  we 
used  a  3  a  median  clip  on  the  data  before  averaging  it  all 
into  a  single  spectrum  with  bins  of  0.05  /mi. 

3.  FINE  STRUCTURE  EMISSION  LINES 

Table  1  lists  the  seven  fine  structure  lines  detected  along 
with  1  a  uncertainties.  The  measured  wavelength  of  each 
line  is  consistent  with  a  225  km  s-1  redshift.  The  t’s  listed 
are  optical  depths  at  each  wavelength  assuming  a  dust 
extinction  model  from  Adams,  Lada,  &  Shu  (1988)  and 
rA  oc  2_1  at  wavelengths  longer  than  ~  50  /mi.  We  calibrate 
this  model  using  the  visual  extinction  from  the  study  of 
near-infrared  hydrogen  recombination  lines  by  Satyapal  et 
al.  (1995).  Assuming  a  foreground  screen  they  found  an  Av 
that  ranged  from  2-12  mag,  with  its  highest  values  near  the 
nucleus  where  Av  ~  10.  We  adopt  this  central  value  of  10 
for  our  visual  extinction,  noting  that  at  the  wavelengths  we 
are  considering  small  errors  in  Av  have  little  effect  on  the 
final  line  fluxes.  The  corrected  flux  is  the  measured  flux 
multiplied  by  exp(r),  appropriate  for  an  external  dust 
screen.  Figure  1  displays  the  entire  spectrum. 

The  fine  structure  lines  we  present  have  been  measured  in 
the  past  by  the  Kuiper  Airborne  Observatory  (KAO ;  Lord 
et  al.  1996;  Petuchowski  et  al.  1994;  Stacey  et  al.  1991; 
Duffy  et  al.  1987;  Lugten  et  al.  1986;  Watson  et  al.  1984), 
but  most  were  taken  on  separate  nights  with  different  cali¬ 
brators.  This  complicates  line  ratios,  because  beam  sizes, 
spectral  resolutions,  and  flux  calibration  schemes  differ. 
Another  advantage  of  this  spectrum  over  KAO  is  the  com¬ 
plete  absence  of  telluric  absorption  features,  which  had 
greatest  effect  on  the  measurements  of  the  [O  i]  63  /mi  and 
[N  n]  122  /mi  lines.  We  have  improved  the  signal-to-noise 
ratio  on  all  lines,  but  the  greatest  increases  were  for  the 
weak  [N  ii]  122,  [N  in]  57,  and  [O  i]  145  jim  lines.  Table  2 
gives  the  KAO  lines  and  their  measured  continuum.  Several 
of  the  KAO  continuum  measurements  provided  are 
published  here  for  the  first  time. 

The  agreement  between  the  KAO  measurements  and  our 
own  fluxes  is  satisfactory,  with  the  most  significant  differ¬ 
ences  in  the  weak  [N  n]  122  and  [O  i]  145  /mi  lines.  Our 
continuum  measurements,  however,  go  from  being  30% 
higher  to  10%  lower  than  the  KAO  values.  Some  of  this 
may  be  due  to  the  smaller  beams  employed  for  KAO  work 
(30"-55"),  but  that  cannot  explain  all  the  discrepancies. 


TABLE  l 


Emission-Line  Fluxes 


Line 

Rest  X 
(/an) 

Flux 

(W  cm-2) 

T 

Corrected  Flux 
(Wcm“j* 

Continuum 

(Jy)b 

[O  hi] . 

51.81 

10.3  ±  0.5  x  10“18 

0.10 

11.3  ±  0.6  x  10“ 18 

1490 

TN  nit . 

57.32 

3.4  ±  0.5  x  10“ 18 

0.09 

3.7  +  0.5  x  10“ 18 

1695 

[O  i] . 

63.18 

17.6  ±  0.5  x  10“ 18 

0.08 

19.1  ±  0.5  x  10“ 18 

1860 

[O  m] . 

88.36 

8.6  ±  0.4  x  10“ 18 

0.06 

9.1  ±  0.4  x  10“ 18 

1970 

[N  ii]  . 

121.90 

1.7  ±  0.3  x  10“ 18 

0.04 

1.8  ±  0.3  x  10“ 18 

1460 

[O  i] . 

145.53 

1.2  ±  0.1  x  10“18 

0.04 

1.2 +  0.1  x  10“ 18 

1155 

[C II] . 

157.74 

13.4  ±  0.1  x  10“ 18 

0.03 

13.8  ±  0.1  x  10“ 18 

1005 

a  Corrected  flux  is  the  measured  flux  corrected  for  extinction. 
b  Continuum  flux  density  has  not  been  corrected  for  extinction. 
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Wavelength(/xrn) 

Fig.  1. — LWS  spectrum  for  M82.  Seven  fine  structure  lines  are  visible, 
as  well  as  an  OH  absorption  at  119  /im.  The  strong  features  seen  around 
43-50  /an  are  unidentified  and  possibly  artifacts  of  the  first  detector. 

The  Infrared  Astronomical  Satellite  (IRAS)  measured  a 
continuum  flux  of  1271  Jy  at  60  gm  and  1351  Jy  at  100  gm 
(Rice  et  al.  1988),  both  lower  measurements  than  those  pre¬ 
sented  here.  Since  the  IRAS  values  represent  the  total  M82 
flux  in  an  area  roughly  10'  x  5'  and  since  the  background 
measured  by  IRAS  and  by  us  in  reference  positions  is  less 
than  2  Jy,  the  reason  for  this  discrepancy  is  unknown. 

While  the  previous  line  work  (Duffy  et  al.  1987 ;  Lugten  et 
al.  1986;  Lord  et  al.  1996)  achieved  calibration  accuracies  to 
the  20%-30%  level,  our  study  provides  the  first  complete, 
medium-resolution,  far-infrared  spectrum,  with  substan¬ 
tially  reduced  uncertainty  (~5%),  no  atmospheric  prob¬ 
lems,  a  standard  aperture,  and  a  single  calibration  scheme. 
Furthermore,  by  using  line  ratios  rather  than  absolute 
fluxes  to  determine  starburst/photodissociation  region 
(PDR)  properties,  most  systematic  uncertainties  should 
cancel. 

4.  MODELING  AND  OTHER  RESULTS 
4.1.  Blackbody  Fitting 

We  fit  our  spectrum  with  a  function  of  the  form  F k  oc 
Bk(T)(  1  —  e_TDusl),  where  we  assume  rDust  oc  2_1.  This  gave 
us  a  good  fit  for  a  temperature  of  48  K  for  the  emitting  dust. 


Log  Wavelength(um) 

Fig.  2. — Blackbody  fit.  The  solid  line  is  actual  data  corrected  for  pre¬ 
dicted  extinction;  the  dotted  line  is  a  best-fit  blackbody  spectrum  with  flux 
density ocBA(T)(l-e_’Du“)  law,  where  t  oc  2_1.  For  comparison,  the  dashed 
line  is  our  best  fit  with  %  oc  X  “  2. 


If  we  then  assume  rDust  =  1  at  ~  2  gm  (from  Av  —  10),  we 
get  a  covering  fraction  of  almost  unity.  Emissivity  laws 
where  r  oc  2-1'5  and  2” 2  were  also  considered.  The  2“ 15 
law  does  not  fit  as  well,  and  the  2” 2  law  showed  large 
deviations.  The  48  K  temperature  is  consistent  with  the 
formula  from  Spinoglio  et  al.  (1995),  which  determines  the 
color  temperature  of  a  galaxy  with  a  2"1  emissivity  law 
from  its  spectral  index  found  using  only  the  flux  densities  at 
60  and  100  gm.  From  our  data  oe6o-ioo  for  M82  is  0.04. 
When  we  substitute  that  into  their  regression  formula, 
Tcolor  =  11.4  x  (a60_ioo  +  4.67)  K,  the  result  is  53  K,  only 
5  K  different  than  a  fit  to  the  entire  far-infrared  spectrum. 

Figure  2  shows  the  fit  with  extrapolation  into  the  sub¬ 
millimeter.  Not  shown  is  Elughes,  Gear,  &  Robson’s  (1994) 
data  point  for  450  gm,  2FA  =  (2.7  ±  0.4)  x  10-17  W  cm-2, 
measuring  the  flux  from  the  central  68"  x  68"  of  M82. 
Inputting  450  gm  into  our  fit  gives  a  2FA  =  5.5  x  10” 17  W 
cm”  2,  a  factor  of  2  high,  indicating,  perhaps,  a  steepening  in 
the  2  dependence  for  the  dust  emissivity  law  in  the  sub¬ 
millimeter.  Klein,  Wielebinski,  &  Morsi  (1988),  fitting  a  few 
points  from  the  infrared,  submillimeter,  and  out  into  the 
millimeter  radio  portion  of  the  spectrum,  did  derive  a 
steeper  dust  emissivity  law,  with  rDust  oc  2-15. 

The  total  flux  for  the  infrared  region  observed  is 
8.15  x  10“ 15  W  cm-2,  which  at  a  distance  of  3.63  Mpc 
corresponds  to  3.2  x  1010  L0.  Assuming  it  continues  short- 


TABLE  2 


KAO  CONTINUA  AND  FLUXES 


KAO  Linea 

KAO  Continuum 

HPBW 

Line 

(W  cm-2) 

(Jy) 

(arcsec) 

Reference 

[O  m]  52  /im . 

9.5  ±  0.7  x  10“ 18 

1207 

48 

1 

[N  m]  57  /im . 

3.9  ±  0.4  x  10“ 18 

1269 

48 

1 

[O  i]  63  fim . 

14.2  +  3.4  x  10“ 18 

1254 

44 

2 

[O  m]  88  /im . 

8.6  ±  0.5  x  10“ 18 

1689 

48 

1 

[Nn]  122  fim . 

2.9+_°0%  x  10“18 

1190 

45 

3 

[O  i]  145  . 

0.84  ±  0.24  x  10“ 18 

1130 

55 

4 

[C  ii]  158  /im . 

14  ±  4.2  x  10“ 18 

1150 

55 

4 

a  All  line  fluxes  have  been  corrected  for  same  extinction  used  in  Table  1. 

References — (1)  Duffy  et  al.  1987;  (2)  Lord  et  al.  1996;  (3)  Petuchowski  et  al.  1994;  (4)  Lugten 
et  al.  1986. 
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ward  of  43  /tm  as  a  blackbody,  the  total  far-infrared  flux 
would  be  9.7  x  10” 15  W  cm-2  or  3.8  x  1010  L0,  in  agree¬ 
ment  with  the  estimate  by  Telesco  &  Harper  (1980).  The 
total  flux  in  the  fine  structure  lines  is  6.0  x  10“ 17  W  cm-2 
or  0.6%  of  the  infrared  flux.  The  [C  ir]/F(FIR)  ratio  in  the 
LWS  beam  is  1.4  x  10“ 3,  similar  to  that  seen  toward  galac¬ 
tic  H  ii  regions  and  other  starburst  and  normal  galaxies 
(Stacey  et  al.  1991). 

4.2.  Emission-Line  Ratio  Modeling 

A  realistic  starburst  of  a  given  age  and  duration  is  inter¬ 
mediate  between  the  simple  cases  of  an  instantaneous  burst 
and  continuous  star  formation.  In  order  to  limit  the  number 
of  free  parameters  and  to  bracket  the  real  situation,  we 
looked  at  two  different  starburst  types,  instantaneous  burst 
and  continuous  star  formation.  For  each  one  our  combined 
starburst  H  n  region  and  photodissociation  region  (PDR) 
model  has  six  free  parameters:  density  in  H  ii  region, 
density  in  PDR,  ionization  parameter  in  H  n  region,  far- 
ultraviolet  flux  in  PDR,  input  starburst  age,  and  input  star- 
burst  upper  mass  limit. 

Table  3  lists  the  measured  line  ratios  of  M82  together 
with  the  predicted  combination  starburst  (H  ii  region)  and 
PDR  model  ratios.  Following  the  methodology  of  Fischer 
et  al.  (1996),  the  H  ii  region  model  was  created  using  the 
program  CLOUDY  by  Gary  Ferland  (Version  90;  Ferland 
1996)  and  input  spectral  energy  distributions  (SEDs)  of  star- 
bursts  described  in  Leitherer  &  Heckman  (1995).  To  test 
different  SED  models,  we  tried  two  upper  cutoff  masses,  30 
and  100  M0,  and  two  burst  types,  instantaneous  bursts  and 
continuous  star  formation,  for  ages  of  1-25  Myr.  All  the 
SED  models  had  solar  metallicity  and  initial  mass  function 
(IMF)  slopes  of  2.35.  The  input  parameters  are  density  and 
the  ionization  parameter,  U,  which  is  defined  as  the  ratio  of 
ionizing  photons  to  hydrogen  atoms  at  the  inner  face  of  the 
cloud : 


where  Q(H )  is  the  number  of  (ionizing  photons)  s_1,  nH  is 
the  total  hydrogen  density  at  the  inner  face  of  the  cloud,  and 
r  is  the  distance  from  the  ionizing  source.  Orion  abundances 


and  dust  grains  were  assumed  for  the  emitting  H  ii  regions. 
The  [O  m]  52  gm/[0  m]  88  /on  ratio  is  sensitive  to  density, 
but  it  is  almost  completely  independent  of  ionization 
parameter  and  model  type.  We  found  the  best  fit  for  a 
density  of  250  cm-3  for  all  models.  The  [N  hi]  57  /un/[N  n] 
122  /on  ratio  acts  in  an  orthogonal  manner  to  our  first 
diagnostic  ratio.  It  is  weakly  dependent  on  density  for  low- 
density  H  ii  regions,  but  it  is  strongly  dependent  on  ioniza¬ 
tion  parameter  and  weakly  dependent  on  model  type.  Here 
we  fit  log  U  varying  from  —  3.1  for  the  30  MG  cutoff  contin¬ 
uous  star  formation  to  —3.5  for  a  3-7  Myr  old  instantane¬ 
ous  burst  with  a  100  Me  cutoff.  These  fits  with  different 
ionization  parameters  demonstrate  Spinoglio  &  Malkan’s 
(1992)  point  that  a  harder  ionizing  spectrum  can  be  some¬ 
what  compensated  for  by  a  modest  decrease  in  U. 

Change  in  input  abundances  made  little  difference  in  our 
derived  starburst  parameters.  We  ran  CLOUDY  models 
with  Orion  abundances,  ISM  abundances,  solar  abun¬ 
dances,  and  twice  solar  metallicity  starburst  abundances, 
where  the  abundances  are  those  defined  in  CLOUDY 
(except  for  the  ISM  abundances,  which  come  from  Sembach 
&  Savage  1996).  No  matter  what  abundances  we  input,  we 
always  derived  basically  the  same  densities  and  ionization 
parameters.  The  quality  of  the  fit  did  vary  somewhat,  but 
this  was  mainly  the  fit  of  oxygen  to  nitrogen  line  ratios 
([O  in]  52  /im/[N  hi]  57  gm,  [O  in]  88  /tm/[N  ii]  122  gm, 
and  so  on),  which  one  would  expect  to  change  when  altering 
the  ratio  of  oxygen  to  nitrogen  abundances.  Orion  abun¬ 
dances  fit  best,  matching  to  within  10%  of  our  measured 
ratios.  ISM  and  solar  abundances  have  only  slightly  differ¬ 
ent  O/N  abundance  ratios,  and  so  they  also  produced  satis¬ 
factory  fits  (within  30%)  of  the  measured  O/N  line  ratios. 
Twice  solar  metallicity  starburst  abundances  have  high 
O/N  abundance  ratios  and  consequently  produced  high 
O/N  line  ratios,  ~  2.5  x  that  observed. 

Of  the  four  model  input  starburst  SEDs  considered,  only 
two  produced  satisfactory  ratio  fits :  the  3-7  Myr  instanta¬ 
neous  burst  with  100  M0  cutoff  and  the  8-25  Myr  contin¬ 
uous  star  formation  with  30  M0  cutoff.  A  fit  was  also 
achieved  with  the  continuous  star  formation  with  100  MQ 
cutoff,  but  only  after  reducing  the  nitrogen  abundance  by  a 
third.  We  can  further  constrain  the  models  by  looking  at 
their  total  bolometric  luminosity.  Using  the  extinction- 


TABLE  3 
Line  Ratios 


Line  Ratio  Measured  Ratio  Instantaneous3  Continuous11 


[O  hi]  52  /an/[N  m]  57/im .  3.0  +  0.4  3.0  3.1 

[O  hi]  52  /im/[0  m]  88  /an .  1.24  ±  0.08  1.26  1.32 

[O  hi]  52  /im/[0  i]  63  /an .  0.59  ±  0.04  0.59  0.59 

[O  hi]  52  /an/[N  n]  122  /an .  6.3  ±  0.5  6.3  6.6 

[O  hi]  52  /a n/[0  i]  145  /mi .  9.4  ±  0.9  9.2  10.5 

[O  hi]  52  /an/[C  n]  158  /an .  0.82  ±  0.04  0.82c  0.82c 

[N  hi]  57  /an/[N  n]  122  /an .  2.1  ±  0.3  2.1  2.1 

[O  i]  63  /a n/[0  i]  145  /an .  15.9  ±  1.4  15.5  17.9 

[O  i]  63  /a n/[C  n]  158  /an  .  1.38  ±  0.03  1.39  1.41 

[N  ii]  122  /an/[N  n]205  /an .  4.1  +  1.4d  4.4  4.9 


3  Instantaneous  burst  model,  3-5  Myr  old  with  100  M0  cutoff  plus  PDR  contribution,  the 
preferred  model. 

b  Continuous  star  formation  model,  8-25  Myr  old  with  30  M0  cutoff  plus  PDR  contribu¬ 
tion.  This  model  was  rejected,  because  it  produced  too  much  luminosity. 

c  This  ratio  is  forced  to  match  exactly  by  the  fitting  method. 

d  ISO  did  not  observe  the  [N  ii]  205  /an  line;  we  use  the  KAO  observation  from  Petu- 
chowski  et  al.  (1994). 
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corrected  Bry  flux  contained  within  an  ISO  LWS  beam 
(from  Satyapal  et  al.  1997)  and  assuming  an  electron  density 
~ 200  cm-3  and  an  electron  temperature  ~ 5000  K,  we 
found  an  ionizing  photon  rate  of  8  x  1053  s_1.  This  is 
similar  to  values  previously  found.  For  instance,  McLeod  at 
al.  (1993)  found  a  rate  of  1.05  x  1054  s_1  using  a  smaller 
beam  (30")  but  different  assumptions  about  extinction.  As 
one  experiment,  we  scaled  the  Leitherer  &  Heckman  (1995) 
models  to  this  ionizing  photon  rate  and  found  that  the 
continuous  star  formation  produces  ~  1 0 1 1  L0,  more  than 
twice  the  measured  far-infrared  luminosity  ( ~  4  x  1010  LQ). 
The  older  instantaneous  bursts  also  produce  too  much 
luminosity,  but  the  more  recent  bursts,  3-5  Myr  old,  only 
produce  (3-5)  x  1010  L0,  more  consistent  with  the 
observed  total  far-infrared  luminosity.  This  makes  the  3-5 
Myr  instantaneous  burst  our  preferred  input  SED.  Small 
errors  in  ionizing  photon  rate  do  not  affect  this  conclusion. 

A  similar  check  was  made  by  comparing  the  Leitherer  & 
Heckman  (1995)  starburst  masses  expected  to  the  dynami¬ 
cal  mass  derived  from  the  mass  model  by  Gotz  et  al.  (1990), 
which  gives  us  1.6  x  109  M0  within  our  ISO  beam. 
McLeod  et  al.  (1993)  also  used  the  Gotz  et  al.  (1990)  mass 
model  but  for  their  smaller  beam  size.  After  they  subtracted 
estimates  for  molecular  mass  and  an  older  stellar  popu¬ 
lation,  they  arrived  at  an  estimate  of  2.5  x  10 8  Me  for  the 
starburst  population.  Rieke  et  al.  (1993)  used  this  mass  limit 
along  with  modeling  to  predict  that  a  lower  mass  cutoff 
must  exist.  Our  preferred  instantaneous  burst  model  pre¬ 
dicts  starburst  masses  of  0.5-1. 3  x  108  Me,  depending  on 
lower  mass  cutoff  (1.0-0.1  M0),  which  does  not  come  close 
to  exceeding  the  dynamical  mass  contained  within  the  ISO 
beam.  Even  if  we  take  the  2.5  x  108  M0  used  by  Rieke  et  al. 
(1993)  as  our  total  allowed  starburst  mass,  we  have  no 
trouble  fitting  the  mass  and  see  no  requirement  for  a  lower 
mass  cutoff,  although  one  might  possibly  exist. 

While  the  CLOUDY  models  fit  the  starburst  lines  well, 
they  produce  little  [O  i]  and  [C  ii]  line  emission.  The 
amount  of  [C  n]  158  fan  line  flux  produced  by  the 
CLOUDY  models  depends  on  the  SED  we  input,  but  it  is 
~  24%-3 1  %  of  the  total  observed  line  flux  for  the  instanta¬ 
neous  burst.  Both  [C  n]  and  [O  i]  lines  are  strongly  pro¬ 
duced  in  PDRs,  which  one  would  expect  to  find  in  the 
interface  between  H  n  regions  and  molecular  clouds.  We  use 
the  PDR  models  of  Kaufman  et  al.  (1998).  These  models 
have  two  free  parameters:  density  and  far-ultraviolet  flux, 
G0,  which  is  expressed  in  units  of  local  Galactic  far- 
ultraviolet  flux,  or  1.6  x  10“  3  ergs  s_1  cm-2.  Once  again 
the  line  ratios  provide  excellent  diagnostics  for  determining 
these  parameters.  The  ratio  of  [O  i]  63  /im/[0  i]  145  /im 
acts  orthogonally  to  the  [O  i]  63  /zm/[C  n]  158  /im  ratio 
over  the  densities  and  far-ultraviolet  fluxes  considered 
(Wolfire,  Tielens,  &  Hollenbach  1990).  If  we  assume  that  the 
remaining  69%-76%  of  the  [C  ii]  not  accounted  for  by  the 
starburst  H  n  region  model  is  coming  from  the  PDRs,  then 
we  have  our  PDR  [O  i]  63  /on/[C  n]  158  /an  ratio  of 
~  1. 6-1.8.  Combining  that  with  the  observed  [O  i]  63 
jum/[0  i]  145  /on  ratio  (~  16)  gives  us  a  density  of  1033 
cm-3  and  a  G0  around  102'8.  These  combined  H  ii  region 
and  PDR  models  fit  not  only  the  line  ratios,  but,  scaled  to 
the  distance  and  size  of  M82,  they  also  fit  the  lines  and 
continuum  seen. 

We  can  check  our  derived  G0  and  PDR  density  using 
another  PDR  diagnostic,  the  ratio  of  the  PDR  line  fluxes  to 
the  infrared  continuum  flux.  For  our  parameters  the 


Kaufman  et  al.  (1998)  models  give  Flux  (C  ii  158  jum  +  C  n 
63  /on  +  O  i  145  /im)/FIR  =  3.9  x  10  “3,  roughly  equal  to 
the  2.9  x  10“ 3  we  get  after  subtracting  our  H  ii  region  line 
contribution.  It  should  be  noted  that  the  Kaufman  et  al. 
(1998)  PDR  models  use  Galactic  ISM  abundances,  not  the 
Orion  abundances  input  into  our  CLOUDY  H  n  region 
models.  However,  for  our  PDR  density  and  G0,  the  differ¬ 
ence  in  using  two  different  abundances  is  negligible. 

Since  the  regions  we  are  modeling  are  in  physical  contact 
with  one  another,  we  should  be  able  to  make  some  contin¬ 
uity  checks  across  the  model  boundaries.  Following  the 
treatment  described  by  Satyapal  et  al.  (1998),  we  start  by 
checking  that  there  is  pressure  equilibrium  between  the  H  ii 
region  and  the  PDR.  The  edge  of  the  H  n  region  has  a 
temperature  of  4000  K,  which  gives  a  P/k  ~  106  cm-3  K. 
Inputing  our  PDR  surface  temperature  of  ~  250  K,  derived 
from  Kaufman  et  al.  (1998),  and  density  gives  a  P/k  ~ 
5  x  105  cm-3  K  or  roughly  the  same  as  our  H  ii  region.  The 
second  boundary  condition  is  that  the  far-ultraviolet  flux 
leaving  the  H  n  region  should  equal  that  entering  our  PDR 
model.  We  measure  the  far-ultraviolet  flux  leaving  the  H  n 
region  by  taking  the  CLOUDY  output  continuum  and  inte¬ 
grating  the  flux  from  6  to  13  eV.  This  depends  strongly  on 
input  SED,  giving  us  G0  =  1 02'2— 102  4  for  the  instantane¬ 
ous  burst  models.  Considering  the  uncertainties,  this  is  con¬ 
sistent  with  our  model  G0  =  102'8. 

Lester  et  al.  (1987)  noted  that  the  [O  hi]  52  /im/[N  hi]  57 
fim  ratio  will  equal  the  actual  0++/N++  ratio  within 
±  50%,  thanks  to  the  similar  critical  densities  for  thermali- 
zation  of  their  respective  emitting  levels  and  ionization 
potentials,  which  are  within  25%  of  each  other.  With  know¬ 
ledge  of  the  electron  density,  which  one  can  get  accurately 
from  the  [O  iii]  52  /im/[0  m]  88  fim  ratio,  a  very  precise 
0  +  +/N++  ratio  can  be  found  that  holds  true  for  a  broad 
range  of  electron  temperatures.  Following  the  calculations 
of  Lester  et  al.  (1987)  and  taking  our  density  of  250  cm-3, 
we  find  a  O4  +/N+  +  ratio  of  3.8.  This  is  slightly  higher  than 
that  reported  by  Duffy  et  al.  (1987)  for  M82,  0  +  +/ 
N+  +  =  3.1.  Both  ratios  are  higher  than  the  average  value  of 
H  ii  regions  toward  the  Galactic  center,  where  <0  +  +/ 
N+  +  >  ~  0.9-1. 4  (Dinerstein  et  al.  1984).  One  possible  con¬ 
clusion  that  could  be  drawn  is  that  the  O/N  ratio  is  truly 
higher  in  M82.  However,  it  is  more  likely  that  this  is  an 
effect  of  the  local  ionizing  conditions.  In  low-ionization 
conditions,  i.e.,  cool  stars,  more  nitrogen  is  doubly  ionized 
than  oxygen,  because  the  ionization  potential  of  N+  is 
smaller  than  that  of  0  +  .  The  difference  in  0  +  +  /N+  +  ratio 
in  the  Galaxy  would  result  from  less  massive  (cooler)  stars 
creating  the  H  ii  regions.  This  is  further  supported  by  our 
model  fit,  where  we  fit  the  0  +  +  and  N++  lines  for  M82 
with  a  high-ionization  starburst  SED  and  Orion  abun¬ 
dances  without  any  oxygen  overabundance.  In  fact,  as  this 
section  mentioned  earlier,  the  more  oxygen-rich  abundance 
models  produced  significantly  worse  fits  to  our  line  ratios. 

Examining  our  best  fit  we  come  within  10%  of  all  the 
ratios  and  closer  on  most.  The  derived  parameters — H  n 
region  density,  PDR  density,  ionization  parameter,  far- 
ultraviolet  flux,  burst  upper  mass  cutoff,  and  burst  age — are 
well  determined,  with  even  small  changes  in  any  of  these 
inputs  creating  ratios  outside  the  uncertainties  of  the  ratios 
measured.  The  majority  of  error  still  present  in  these  fits 
comes  from  the  models  themselves  and  the  assumptions  we 
made,  not  from  the  inaccuracy  in  our  line  ratio  measure¬ 
ments. 
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We  can  compare  this  work  with  previous  models.  Spino- 
glio  &  Malkan  (1992)  also  fit  far-infrared  lines  with  a  com¬ 
bination  H  ii  region  and  PDR  model.  While  achieving 
similar  H  n  region  densities,  their  M82  model  had  a  much 
larger  ionization  parameter  (10-25).  This  is  entirely  due  to 
the  hardness  of  the  input  SED.  The  Leitherer  &  Heckman 
(1995)  SEDs  are  significantly  harder  than  the  SED  Spino- 
glio  &  Malkan  (1992)  used.  The  Spinoglio  &  Malkan  (1992) 
PDR  parameters  were  also  different  (G0  =  104,  n  =  104-5 
cm-3),  but  this  can  be  mainly  attributed  to  the  newer  PDR 
models  of  Kaufman  et  al.  (1998).  The  Lord  et  al.  (1996)  PDR 
parameters  for  M82  (G0  =  103,  n  =  104  cm-3)  come  closer 
to  our  derived  values,  but  they  also  differ  because  of  the 
newer  models  used. 

Satyapal  et  al.  (1997)  and  Reike  et  al.  (1993)  previously 
created  models  of  the  M82  starburst,  both  finding  typical 
starburst  ages  of  around  107  yr.  Reike  et  al.  (1993)  com¬ 
pared  many  observational  parameters  ( MK ,  mass,  lumi¬ 
nosity,  CO  index,  UV  photons),  fitting  M82  with  two 
separate  starbursts,  the  first  being  13-30  Myr  ago.  Satyapal 
et  al.  (1997)  looked  at  individual  starburst  clusters,  measur¬ 
ing  their  CO  indices  and  Bry  equivalent  widths.  Modeling 
those  features  for  each  individual  cluster,  they  found  a  range 
in  age  of  6  x  106  yr,  implying  an  outward  propagating  star- 
burst.  For  the  sake  of  simplicity,  we  did  not  attempt  to 
combine  multiple  models  with  different  starburst  ages,  but 
we  note  many  of  the  individual  starburst  cluster  ages  found 
by  Satyapal  et  al.  (1997)  were  around  5  Myr. 

One  unknown  that  could  throw  off  the  models  is  the 
degree  of  [O  i]  63  /un  self-absorption  in  M82.  Stacey  et  al. 
(1983)  first  observed  the  145  /un  [O  i]  line  from  the  PDR 
associated  with  the  Orion  A  H  n  region  and  discovered  that 
the  [O  i]  63,  145,  and  [C  n]  158  /im  line  intensity  ratios 
indicated  optically  thick  (t  ~  2)  emission  in  the  63  /an  line. 
Tielens  &  Hollenbach  (1985)  models  also  indicated  that  this 
should  be  so.  Later,  through  more  sophisticated  modeling 
involving  the  [13C  ii]  line  and  more  complete  mapping  in 
the  [C  ii]  and  [O  i]  lines,  Stacey  et  al.  (1993)  estimated  that 
the  [O  i]  63  /im  has  an  optical  depth  ~  3  in  the  Orion  PDR. 
The  [O  i]  63  /an  line  is  seen  in  absorption  in  Arp  220 
(Fischer  et  al.  1998)  and  in  Galactic  star  formation  regions 
(Poglitsch  et  al.  1996;  Baluteau  et  al.  1997).  In  addition, 
Hermann  et  al.  (1997)  explained  an  unusually  low  [O  i]  63 
/im/[0  i]  145  /an  ratio  in  the  Dark  Lane  region  of  the 
Orion  Molecular  Cloud  as  a  likely  site  of  [O  i]  63  /an 
self-absorption. 

One  would  expect  [O  i]  63  /an  absorption  to  occur  where 
cold  foreground  oxygen  lies  between  us  and  an  emitting 
source.  However,  the  gas  density  must  exceed  the  critical 
density  for  collisional  de-excitation  (ncrit  ~  5  x  105  cm-3) 
in  order  for  the  absorbed  photons  to  be  removed  from  our 
beam.  At  lower  densities  the  photons  are  mainly  scattered, 
with  little  reduction  in  the  observed  line  luminosity.  In  our 
Galaxy,  such  dense  clouds  rarely  cover  more  than  a  small 
fraction  of  the  [O  i]-emitting  region.  Since  the  ISM  in  most 
starburst  and  normal  galaxies  seems  to  be  constructed  of  a 
superposition  of  many  such  clouds,  which  usually  do  not 
cover  each  other  simultaneously  in  both  velocity  and  area, 
the  small  regions  of  [O  i]  self  absorption  should  not  be 
important.  A  notable  exception  is  Arp  220  (Fischer  et  al. 
1998),  where  it  is  likely  the  FIR  source  is  largely  enveloped 
by  cold,  dense  clouds  with  abundant  neutral  oxygen.  There 
is  no  reason  to  think  of  M82  as  such  an  extreme  case,  and  in 
fact  our  analysis  of  high  signal-to-noise  ratio  Fabry-Perot 


scans  of  the  [O  i]  63  /tm  line  does  not  indicate  any  self¬ 
absorption.  Any  significant  [O  i]  self-absorption  would 
alter  the  PDR  modeling  by  increasing  the  PDR  density. 

There  is  one  more  important  forbidden  line  in  the  far- 
infrared,  the  [N  ii]  205  /tm  line,  which  lies  just  beyond  the 
LWS  limit  of  196.7  microns.  Petuchowski  et  al.  (1994)  mea¬ 
sured  it  using  the  KAO,  and  our  model  does  predict  its 
relative  flux.  We  include  both  on  our  line  ratio  table  and 
note  they  match  to  within  the  uncertainty  of  the  KAO  flux 
measurements.  The  [N  ii]  line  ratio  is  the  best  density  tracer 
for  low-density  ( n  <  500  cm-3)  H  ii  regions,  so  the  accurate 
reproduction  of  this  line  ratio  lends  credence  to  our  esti¬ 
mate:  ne  ~  250  cm-3. 

4.3.  Cross  Scans  of  the  C  ii  and  O I  Lines 

In  addition  to  the  full  scan  made  at  the  center  of  M82,  we 
collected  cross-scans  of  the  [C  ii]  158  and  [O  i]  63  /an  lines 
using  the  LWS02  Astronomical  Observation  Template 
(AOT).  The  maps  consist  of  two  cuts  across  M82,  one 
across  the  major  axis  with  a  P.A.  =  55°  and  the  other  across 
the  minor  axis,  perpendicular  to  the  first.  Seven  obser¬ 
vations  were  made  in  each  cut,  starting  150"  from  M82’s 
center  and  spaced  50"  apart.  Each  line  represents  an  inte¬ 
gration  of  3  s  (spectral  element)  - 1.  The  two  plots  in  Figure  3 
show  the  relative  strengths  of  the  emission  lines  along  the 
two  axes.  Only  upper  limits  were  found  for  the  [O  i]  63  /tm 
line  at  the  points  150"  from  M82’s  center  and  are  therefore 
not  included  on  the  map.  The  major  axis  cut  runs  from 
southwest  to  northeast,  and  the  minor  axis  cut  runs  from 
northwest  to  southeast. 

Both  lines  weaken  away  from  the  center,  but  while  the 
[C  ii]  158  /tm  line  fall-off  is  symmetrical  about  the  minor 
axis,  it  is  not  symmetrical  about  the  major  axis;  the  [O  i]  63 
/tm  line  is  not  symmetrical  around  either  axis.  This  asym¬ 
metry  is  not  that  surprising.  Previous  maps,  such  as  those  in 
Ha  (Waller,  Gurwell,  &  Tamura  1992),  the  near-infrared 
(Satyapal  et  al.  1995),  and  the  submillimeter  (Hughes  et  al. 
1994),  show  an  excess  in  flux  to  the  west.  Stacey  et  al.  (1991) 
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Fig.  3. — Plots  of  the  strengths  of  two  lines  across  the  major  (top)  and 
minor  ( bottom )  axis  of  M82.  The  solid  line  is  the  [C  n]  158  pm  line,  and  the 
dashed  line  is  the  [O  i]  63  pm  line.  The  LWS  beam  is  represented  by  the 
dotted  line. 
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Fig.  4. — Plot  of  [O  i]/[C  n]  line  ratio  across  M82.  The  solid  line  is  the 
ratio  across  the  major  axis,  and  the  dotted  line  is  the  ratio  along  the  minor 
axis. 


did  similar  cuts  in  the  [C  n]  158  pm  line,  which  also  showed 
this  slight  asymmetry  in  the  major  axis.  Lord  et  al.  (1996) 
did  detailed  modeling  of  [O  i]  63  pm  line  profile,  attributing 
the  majority  of  the  line  flux  to  two  off-center  hot  spots,  with 
the  western  spot  the  brighter  of  the  two.  We  examined  the 
far-infrared  continuum  in  our  cross  scans  and  also  found 
more  flux  to  the  west. 

We  measured  the  spectral  index  a  between  60  and  100 
pm,  where  Fv  oc  va,  at  every  spot  in  our  scans.  This  a6o-ioo 
is  largest  50"  to  the  west,  where  it  is  ~0.4,  not  at  M82’s 
center,  where  a6o-too  's  ~  0. 1 .  This  higher  a60_100  implies  a 
higher  average  temperature,  further  indication  of  a  strong 
hot  spot  to  the  west.  In  addition,  we  also  looked  at  the 
spectral  index  a  between  100  and  175  pm,  where  we  were 
interested  in  testing  whether  our  single  temperature  black- 
body  with  emissivity  oci  1  would  fit  further  outward  from 
the  galaxy’s  center.  The  a100_175  at  the  hot  spot  to  the  west 
is  ~  1.6,  which  is  consistent  with  our  blackbody  fit. 
However,  the  a,  00-175  measured  to  the  east  and  farther 
outward  in  the  galaxy  is  smaller  than  expected,  possibly 
indicating  a  growing  cold  dust  component. 

Most  interesting  is  the  variation  in  the  [O  1]  63  pm/[C  11] 
158  pm  line  ratio,  visible  in  Figure  4.  Assuming  the  two  lines 
come  mostly  from  PDRs,  tracing  the  [0 1]  63  yum/[C  11]  158 
pm  line  ratios  will  tell  us  something  about  conditions  in 
M82  on  a  large  scale.  The  ratio  shows  a  maximum  roughly 
50"  off  center  to  the  northwest  and  southwest  of  the  galaxy’s 
center  and  shows  a  steady  decrease  outward.  According  to 
the  models  of  Kaufman  et  al.  (1998),  this  ratio  decreases 
with  PDR  density  and  far-ultraviolet  field  strength.  The 
general  trend  of  decreasing  [O  1]  63  pm/[C  11]  158  /on  ratio 
with  distance  from  M82’s  center  is  then  understandable. 
One  would  expect  both  density  and  far-ultraviolet  field  to 
decrease  with  distance  from  the  galaxy’s  center.  The  off- 
center  maximum  in  the  [O  1]  63  pm/[C  11]  158  /on  ratio 
must  indicate  where  conditions  of  high  density,  high  G0,  or 
both  exist. 

The  [C  11]  158  /on  line  flux  also  appears  more  extended 
than  the  [O  1]  63  /on  line  flux,  possibly  indicating  its  associ¬ 
ation  with  the  filamentary  features  observed  in  optical 
recombination  lines,  as  previously  suggested  by  Stacey  et  al. 
(1991).  Taking  the  formula  for  [C  11]  intensity  originating 


from  ionized  regions  from  Madden  et  al.  (1993),  we  get 


7c+(ionized) 


fcv  _A _ gjg, _ 

4?r  ncrit  1  +  Ugjg,)  +  l](ne/ncrit) 


Xc+  EM  , 


(2) 


where  A  is  the  spontaneous  emission  coefficient  for  the 
2P3/2-2Pi/2  transition  (2.36  x  10-6  s_1),  gjgt  is  the  ratio  of 
statistical  weights  in  the  upper  and  lower  levels  (2),  ncrit  is 
the  critical  density  (~35  cm-3),  Xc+  is  the  abundance  of 
C+  relative  to  hydrogen  ( ~  3.3  x  10-4),  and  EM  is  the 
emission  measure.  Assuming  electron  density  is  low  com¬ 
pared  to  the  critical  density,  we  only  need  a  measurement  of 
the  emission  measure  to  derive  the  expected  flux  from  the 
extended,  diffuse  ionized  gas  of  M82.  We  derive  the  EM 
from  the  6  cm  maps  of  Seaquist  &  Odegard  (1991), 
assuming  the  majority  of  the  emission  comes  from  thermal 
free-free  emission.  At  roughly  50"  off  the  minor  axis  we  find 
an  EM  ~  1300  cm-6  pc,  giving  us  an  expected  intensity  for 
C+ofl.7  x  10” 4  ergs  s” 1  cm-2  sr_1  or  an  expected  flux  in 
our  ISO  LWS  beam  of  2.0  x  10-18  W  cm-2.  This  number 
is  consistent  with  that  observed,  indicating  that  the  major¬ 
ity  of  the  extended  C+  emission  is  coming  from  the  diffuse 
ionized  medium. 


4.4.  OH  absorption 

The  only  certain  absorption  feature  observed  in  the  spec¬ 
trum  of  M82  is  due  to  the  two  lambda-doubling  rotational 
transitions  of  OH  at  ~  119.4  pm,  which  are  unresolved  at 
the  LWS  grating-mode  resolution.  These  are  the  transitions 
from  the  OH  ground  state  of  2n3/2  J  =  3/2  to  the  next 
highest  energy  state,  27t3/2  J  =  5/2,  which  have  been  seen 
before  in  emission  in  the  KL  nebula  (Storey,  Watson,  & 
Townes  1981)  and  in  absorption  toward  the  Galactic  center 
(Genzel  et  al.  1985).  No  other  OH  lines  are  observed  in 
M82.  There  is  a  dip  at  the  location  of  the  53  pm  transitions, 
but  it  is  the  same  scale  as  the  noise  and  must  be  treated  as 
an  upper  limit.  Table  4  lists  the  1  a  upper  limits  for  the  next 
most  likely  OH  lines. 

The  equivalent  width  of  an  optically  thin  line  is 


W  = 


gj  87 ic 


(3) 


where  Atj  is  the  Einstein  coefficient  for  spontaneous  emis¬ 
sion  between  two  levels,  i  and  j;g  is  their  statistical  weights; 
X0  is  the  transition’s  line  center;  and  N  is  the  column  density 
of  material  (Aller  1960).  Taking  the  A  coefficients  from  Des- 
tombes  et  al.  (1977),  one  can  calculate  column  densities  and 
column  density  upper  limits  for  the  lines.  Three  of  the  tran¬ 
sitions,  53,  79,  and  119  pm,  are  all  from  the  ground  state, 
but  we  only  detect  the  119  pm  line  because  the  A- 
coefficients  for  53  and  79  pm  give  equivalent  widths  20-80 
times  less  than  that  of  the  119  pm  transition,  assuming  the 


TABLE  4 
OH  Lines 


Line  Equivalent  Width  Atj 

(;im)  Transition  (/im  )  (s-1) 


53 .  2ji3/2  T  =  3/2  ->  7t1/2  2  J  =  3/2  <0.003  0.04 

79 .  2ji3/2  J  =  3/2  -►  tc1/2  2  J  =  1/2  <  0.004  0.033 

84 .  2n3/2J  =  5/2->k3I22  J  =  7/2  <0.004  0.5 

119 .  27r3/2  7  =  3/2-*  7t3/2  2  J  =  5/2  0.022  +  0.004  0.13 

163 .  2tc1/2  J  =  1/2  tc3/2  2  J  =  3/2  <  0.005  0.054 
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transitions  are  optically  thin.  This  is  consistent  with  the 
upper  limits  observed.  The  lack  of  the  84  and  163  /tm  lines, 
despite  having  higher  cross  sections,  demonstrates  that  the 
absorbing  OH  gas  is  almost  entirely  in  its  ground  state.  The 
column  density  of  OH  in  its  ground  rotational  state  is 
~4.2  x  1014  cm-2.  Further  analysis  will  be  presented  in  a 
later  paper. 

Besides  the  119  /an  OH  absorption  lines  and  the  seven 
emission  lines,  we  identify  no  other  lines  or  other  features  in 
the  spectrum.  We  should,  however,  point  out  some  visible 
features  in  the  spectrum  that  may  or  may  not  be  real.  The 
large  bumps  at  the  shortest  wavelengths,  43-50  n m,  are 
present  in  other  LWS  spectra  (Arp  220;  Fischer  et  al.  1998). 
These  features  appear  to  be  an  artifact  of  that  detector, 
SW1,  but  they  are  being  investigated  further.  Two  other 
apparent  features  at  56  and  112  /tm,  result  from  unexpected 
HD  lines  in  the  Uranus  calibration  spectrum.  The  56  /an 
line  is  particularly  troublesome,  because  it  confuses  the  situ¬ 
ation  on  the  short  wavelength  side  of  the  [N  m]  57  /an  line, 
where  there  is  a  hint  of  an  absorption,  possibly  a  p-H20 
line.  Harvey  et  al  (1998)  have  examined  the  7S0/LWS  M82 
spectrum  for  evidence  of  atomic  hydrogen  recombination 
lines  in  this  range.  None  were  detected.  The  upper  limits, 
constrained  by  their  very  low  line-to-continuum  ratio,  lie 
within  the  range  of  expected  values  set  by  the  known  centi¬ 
meter  and  decimeter  lines  and  a  spontaneous  emission 
process.  Other  potential  features  not  already  mentioned  are 
absorptions  at  65.1,  66.4,  and  149.2  /tm.  Further  analysis 
and  possible  identification  will  come  in  a  later  paper.  Any 
undetected  lines  present  in  the  M82  spectrum  would  have 
to  have  fluxes  less  than  ~5  x  10-19  W  cm-2  at  50  /tm  to 
~2  x  10” 19  W  cm-2  at  180  /um. 

5.  CONCLUSIONS 

We  have  obtained  an  LWS  full  grating  scan  of  the  pro¬ 
totypical  starburst  galaxy,  M82.  This  spectrum  reveals 
seven  bright  fine  structure  lines,  previously  detected  from 
the  KAO  but  here  detected  at  the  same  time,  with  the  same 
instrument,  with  roughly  the  same  beam,  and  without  any 
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obstructing  atmosphere.  Therefore,  the  relative  calibration 
of  these  lines  is  superior  to  that  used  for  prior  analysis.  Also 
reported  is  the  discovery  of  an  absorption  line  of  OH  from 
its  lowest  state.  The  upper  limits  for  other  possible  OH  lines 
are  enough  to  determine  that  most  of  the  absorbing  OH  is 
in  its  ground  state.  The  whole  spectrum  is  well  fitted  by  a 
single  temperature  dust  component  with  an  optical  depth 
^Dust ,QC  giving  a  total  infrared  flux  of  3.8  x  1010  L0. 

We  fitted  the  infrared  line  ratios  of  M82  with  a  six- 
parameter  combination  H  ii  region  and  PDR  model.  The 
parameters  for  the  best  fitting  model  are  H  ii  region 
density  =  250  cm-3,  U  =  10~3S,  PDR  density  =  1033 
cm-3,  G0  =  102  8,  upper  mass  cutoff  =  100  M0  and 
age  =  3-5  Myr.  The  [C  ii]  158  jum  line  comes  ~25%  from 
the  H  ii  regions  and  ~75%  from  PDRs.  In  our  model,  the 
starburst  mass  is  ~  (0.5-1. 3)  x  10 8  M0,  depending  on  the 
lower  mass  cutoff  chosen.  This  is  a  significant  percentage  of 
the  molecular  mass  (~2  x  108  Me,  Wild  et  al.  1992),  indi¬ 
cating  the  M82  starburst  will  not  be  able  to  continue 
producing  stars,  through  another  instantaneous  burst  or 
otherwise,  for  much  longer  before  running  out  of  mass.  The 
LWS  cross-scan  data  support  the  model  of  strong,  off- 
center  hot  spots  seen  in  the  mid-infrared  (Telesco  et  al. 
1991)  and  modeled  from  velocity  profiles  by  Lord  et  al. 
(1996).  Our  fit  for  a  3-5  Myr  instantaneous  starburst  is  a 
simplification  of  the  true  situation,  which  is  several  regions 
with  different  starburst  ages.  The  single  burst  we  modeled 
should  be  dominated  by  the  brightest  and  most  recent 
bursts,  which  are  presumably  the  hot  spots  themselves. 
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